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ABSTRACT 

The presence of a circumnuclear stellar disk around Sgr A* and megamaser systems near other black 
holes indicates that dense neutral disks can be found in galactic nuclei. We show that depending on 
their inclination angle, optical depth, and spin temperature, these disks could be observed spectro¬ 
scopically through 21 cm absorption. Related spectroscopic observations of Sgr A* can determine its 
HI disk parameters and the possible presence of gaps in the disk. Clumps of dense gas similar to the 
G2 could could also be detected in 21 cm absorption against Sgr A* radio emission. 


1. INTRODUCTION 


The presence of a disk of massive young stars ar ound 
the Galactic Centre (Levin and Beloborodo v |2003|) and 
various megamaser systems (iMiyosETeF^i.ll 19951 [Moran 


et al. 2007 


Herrnstein et al.]|2005| Kuo et al. 2bll[ ) in¬ 


dicates that dense neutral disks can be found in galactic 
nuclei. These disks could be produced by the tidal dis¬ 
ruption of molecular clo uds passing close to the centra l 
supermassive black hole ( Yusef-Zadeh and Wardle|2012 ). 
The disk could be an outer extension of the hot accretion 
disk closer to the black hole, and possess a large amount 
of neutral hydrogen. Indeed, a significant portion of the 
H 2 O megamasers show X-ray abso rption with large col¬ 
umn densities of 10^^ — 10^® c m~^(I^eb||2008 ). Fu rther, 
the discovery of the G 2 cloud (Cillessen et al.|2(J14 ) indi¬ 
cates that dense clouds of self-shielded neutral hydrogen 
may exist around Sgr A*. 

Here we present a method to observe HI disks in the in¬ 
termediate region of 10^ —10"^ Schwarzschild radii around 
nuclear black holes through their 21 cm absorption. The 
inherent brightness of such a disk is too weak to be de¬ 
tectable, however the black hole’s inner accre tion flow 
emits synchrotron radiation (Bower et al. 2014) that acts 
as a background source upon which the neutral systems 
can be seen in absorption. We examine how observations 
of the absorption spectrum allows one to determine the 
parameters of the HI absorber without the need to spa¬ 
tially resolve the system. 


2. 21 CM ABSORPTION PROFILE 

The accretion flow near Sgr A* is optically-thick to 
synchrotron self-absorption at a wavelength of 21 cm, 
and this photosphere illuminates the outer HI disk. At 
a distance D, the photosphere intensity yields a flux 
of 
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where lyQ is the photospheric specific intensity, (r, 6) are 
coordinates on the disk plane centered on Sgr A*, and the 
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Fig. 1.— A sketch of the geometry from the side (top image) 
and as seen by the observer (bottom image). The black hole is 
surrounded by an optically thick photosphere with a radius 
and an HI disk. The solid arrow in the top image indicates the 
direction towards the observer. The inclination angle i is the angle 
between the observer and the normal to the disk. The dotted 
arrows in the top image depict photons leaving the photosphere in 
the direction of the observer. In their path towards the observer, 
these photons encounter the HI disk and their intensity is reduced 
by absorption. The region of the HI disk that is illuminated by 
these photons is shaded green in both images. The outer radius 
of this region, as seen from the observer’s point of view (bottom 
image), is referred to as Rmax- 


integral is over the area A of the disk that is illuminated 
by Sgr A*’s photosphe re. The freq uency dependent op¬ 
tical depth is given by ( Loeb|[ 2008 1 , 


— sin i cos 9 , 
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where M is the black hole mass, i the disk inclination 
defined as the angle between the observer and the normal 
to the disk (refer to Figure!^, A 21 = 2.85 x 10“^® s“^ 
is the Einstein coefficient of the 21 cm transition, E/k = 
0.068 K is the 21 cm transition energy, Nh is the HI 
column density in the disk, v is the frequency, and T 5 is 
the gas spin temperature. For our calculations, the line 
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function </) is a gaussian, 
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with (T ^ 1/21 = 1-4 X 10® Hz. For a black hole with 
a photosphere radius Rph, the area A can be evaluated 
geometrically. Assuming that the disk is larger than the 
photosphere, the amount of disk that is illuminated by 
the photosphere is the half-circular area spanned by the 
photosphere projected onto the disk. Only half a circle 
is required because the other half the disk is positioned 
behind the photosphere. For a disk with an inclination 
i, the projection gives 
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is the outer radius of the region in the disk that partic¬ 
ipates in absorption (see Figure]^. The absolute value 
of the term proportional to can be written as 
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In general, Njj and T 5 may depend on the radial coor¬ 
dinate r. In addition, 1,^^ also depends on both r and 9 
due to limb darkening. 


2 . 1 . Homogeneous disk with no limh darkening 


For simplicity, we focus our attention to a homoge¬ 
neous disk where N^, Tg, and are constants. From 
equation (§, 
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Besides the multiplicative factor the problem 

possesses only two free parameters: the disk inclination 
angle, i, and the radius of the photosphere, Rph- Al¬ 
though the intrinsic size of of Sgr A* was never measured 
at a wavelen gth as long a.s 21 cm. we extrapolate from 

(20141 and find Rph ~ 10 ^^ cm. 


the results of 


Bower et al. 


The depth of ttie line protiie relative to the continuum. 
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Fig. 2. — The neutral disk as seen by the observer. Due to the 
Doppler shift, an observer sees a line of sight velocity profile that 
traces a dipole on the neutral disk (The black hole is located at 
x,y = 0). The color describes the magnitude of the line of sight 
velocity; areas with larger (lower) redshifts are redder (bluer). The 
value in the color bar is normalized to y/GM sin-i/r. The boundary 
of the disk that is illuminated by the photosphere, Rmaxi is plotted 
in dashed lines for a variety of photosphere radii, Rph (top) and 
a variety of inclinations, i (bottom). The outermost Rmax in the 
top panel corresponds to Rph = 10^^ cm, and the successively 
smaller contours are for Rph = 5 x 10^^ cm and Rph = 3 x 10^^ 
cm. The outermost Rmax in the bottom panel corresponds to 
i = 15°, and the successively smaller contours corresponds to 2 = 
30°, 45°, 60°, Sc 75°. The x and y axes are Cartesian coordinates 
on the disk, in units of 10^^ cm. The portion of the disk that is seen 
in absorption is within 0 < 0 < tv and Rph ^ < Rmax- Increasing 

Rph samples a wider region of the disk, but a larger part of the high 
velocity portion of the disk close to the black hole is hidden behind 
the photosphere. This results in an absorption profile that is deeper 
and thinner. A more edge-on (higher i) disk samples a wider region 
of the disk. Furthermore, the shape of the illuminated portion of 
a more edge-on disk covers higher velocity features, generating a 
deeper and wider absorption profile (see Figure]^. 

SF^/F^g, is shown in Figure]^ for Nh = 10®^ cm“® and 
Tg = 8 X 10^ K. Note that SF^, oc Nn/Tg more gen¬ 
erally. The distance to the Galactic center is taken as 
D = 8 kpc. A physical understanding of the absorption 
profiles can be obtained by looking at the line of sight 
velocity structure of the disk as shown in Figure Due 
to the Doppler shift, an observer would detect a hne of 
sight velocity that varies like a dipole on the disk (con¬ 
tours of ~ cos9/^/r). The portion of the disk that is 
illuminated by the photosphere, i.e. between 0 < 6 * < tt 
and Rph < r < Rmax, is seen in absorption. At every 
point in the illuminated disk (parameterized by the disk 
coordinates r and 9), there is a corresponding drop of 
flux in the frequency profile at i/ = GM sin i cos 9 jrc^ 
due to the absorption at that point. 

The change in the shape of the absorption prohle when 
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Fig. 3. — Absorption profile of the Sgr A* photosphere for 
= 10^^ cm and an inclination angle of 2 = 7r/4 (solid), 
= 2 X 10^^ cm and i = 7r/4 (dashed), as well as Rph = 10^^ 
cm and i = TrjS (dotted). Vr denotes the line-of-sight velocity. 
Since the disk is behind the photosphere when r < Rph^ a larger 
photosphere samples less of the high velocity regions of the disk, 
resulting in a thinner profile. However, a larger photosphere illu¬ 
minates a larger disk area, resulting in more absorption. Increasing 
Rph therefore results in a deeper, but thinner profile. Edge-on ge¬ 
ometry implies a larger portion of the disk is illuminated by the 
photosphere, therefore increasing the angle i results in a deeper 
profile. A disk with higher i samples higher velocity regions, re¬ 
sulting in a wider profile. These extra high velocity regions are 
sampled at the edges of the illuminated disk, in contrast to the 
high velocity regions sampled in systems with larger Rph-) where 
the extra regions are located closer to the black hole (see Figure 

i- 

one changes Rph and i can be explained by the shape of 
RrnaxiRph, i) On the neutral disk. As Figure]^ shows, the 
outer bounds of the illuminated disk increases with in¬ 
creasing Rph- This increases the area of the illuminated 
disk, generating a deeper absorption profile. However, as 
Rph increases, a larger portion of the inner disk is hid¬ 
den behind the photosphere. This removes high velocity 
components from the int^ral, causing a thinner spectral 
profile. Indeed, Figure demonstrates that increasing 
Rph produces a deeper and thinner absorption profile. 
A more edge-on (higher i) disk also possesses a larger 
illuminated disk. In this case, the inner bound is not 
changing, thus the wider illuminated disk also contains 
more high velocity portions. This generates a deeper and 
wider profile, as seen in Figure]^ 


2 .2. Disks with gaps 

Using the formalism of equation Q we can also cal¬ 
culate the resulting 21 cm absorption profile in the case 
where there are gaps in the HI disk. This is done by 
simply breaking the radial integral into parts, 
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where I is the integrand of equation Q, Rm is the in¬ 
ner radius of the gap, and Rout is the gap’s outer radius. 
The integral can be broken to more pieces if more gaps 
are present. An example profile for a disk with a gap is 
shown in Figure]^ A disk with a gap lacks absorption 
on the region where the gap is located. Based on the 
dipolar contours of Figure the velocities correspond- 


Fig. 4.— 21 cm absorption profile for a disk inclined by 80 de¬ 
grees and Rph = 2 x 10^^ cm with a gap of width ~ Rph (solid) 
overplotted against a disk with no gap (dashed). Vr denotes the 
line-of-sight velocity. 


ing to a portion between Rin < r < Rmax is excluded 
from the integral, generating a profile with visible wings. 
The detection of a gap can be used to constraint the 
presence of an intermediate mass black hole companion 
of Sgr A* orbiting with a semimajor axis of ^ 10^ — 10^ 
Schwarzschild radii. 


2.3. Orbiting dense cloud 

Another source of HI is a dense cloud that orbits in 
front of the 21 cm photosphere. This dense spot could 
be a feature in the neutral dis k or a clump of self-sh ielded 
gas similar to the G2 cloud ( Gillessen et al. [2014 |. The 
resulting profile will be like that of Figure]^ but with an 
extra absorption feature located where the spot resides in 
frequency space, as shown schematically in Figure!^ The 
core of this feature travels in i/-space as the cloudorbits 
the black hole. If the spot orbits in a perfect circular 
orbit, the frequency position of the center of the feature, 
Vc, obeys 


z/c(t) = l.dxlO^Hz 


, 1 [GM . . / GM, 

where Tc is the radial position of the dense spot, and t 
the time coordinate relative to an arbitrary initial time 
to- The resulting effect in frequency space is a feature 
that oscillates sinusoidally around the 21 cm rest frame 

frequency. The amplitude of the oscillation is sin i 

and the temporal frequency is Since both of these 

are observables, one can use them to measure simultane¬ 
ously the orbital radius and the black hole mass up to a 
factor of sinf. Plunging orbits can be treated by adding 
a time dependence to Tc, 
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where a is the semimajor axis of the orbit and 

a(l — e^) 
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Fig. 5. — A schematic illustration of an orbiting dense spot. As 
the spot orbits, it covers different portions of the line of sight ve¬ 
locity structure. This resulted in extra absorption on top of the 
neutral disk profile. The feature corresponding to the dense spot 
travels in v space according to in equation § , while the HI 

disk profile is stationary. For a cloud on a circular orbit, the spot 
absorption feature oscillates sinusoidally in frequency space. The 
amplitude and frequency of the oscillation can be used to simulta¬ 
neously measure the spot orbital radius and the black hole mass up 
to a factor of sinz. Note that a dense spot can exist in the absence 
of a disk. 


As such, the detectability of the signal would depend on 
N}i and Ts of the observed system. In particular, for 
i = 80 degrees disk, the signal to noise is given by 
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As we have no constraints on the values of Nh and Tg 
at the Galactic Center, the fiducial values of and 
Tg considered are taken from observations of megamaser 
systems and theoretical considerations in our G alactic 
Center ( Loeb|2008 Yusef-Zadeh and Wardle|2dT2 ). Note 
that the hducial value ot Njg m our calculation is modest, 
and column density in excess of 10 ^^ cm“^ might be pos¬ 
sible. The spin temperature, T 5 , is even less constrained 
than Nii\ while it could be coupled to the actual gas 
temperature, our model considers Tg as a free parameter 
that is not necessarily equal to the gas temperature. Our 
method promises to constrain the combination of the two 
parameters, Nn/Tg. 


3. CONCLUSIONS 


with e being the orbital eccentricity. In addition to the 
inclination angle, the orbit possesses three parameters: 
M, a, and e (or equivalently, the orbital energy, angular 
momentum, and the black hole mass). If the mass of the 
black hole is known, the oscillation amplitude and fre¬ 
quency can be used to determine the orbital parameters 
of the cloud up to the unknown inclination factor. 


2.4. Detectability 

The signal to noise ratio of such a system is given by 
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where SEFD is the System Equivalent Flux Density 
of the telescope, Ai/ the bandwidth, to the observing 
time, and the factor of results from the use of dual 
polarization observations. If Sgr A*’s intrinsic (unab¬ 
sorbed) 21 cm flux density is ~ I Jy, Nh = 10^^ cm“^, 
Tg = 8 X 10^ K, i = 80 degrees, and Rph = 10^® cm, we 
find that SF^, ^ 0.6 mJy. The Square Kilometer Array 
(SKA) possesses a collecting area of 10^° cm^. A system 
temperature of ~ 50 K will then give SEFD = 0.3 Jy. 
Assuming Aiz ~ 10 MHz, we obtain S/N > 10 over 30 
minutes of observations. We also note that the signal to 
noise ratio scales with column density, spin temperature, 
and observation time as 
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We have calculated the 21 cm absorption profile cor¬ 
responding to a disk of neutral hydrogen around black 
holes. Spectroscopic measurements of a disk absorption 
profile can be used to determine the disk parameters, or 
to ascertain the presence of gaps in the disk even when 
the black hole-disk system is spatially unresolved. Fur¬ 
thermore, we demonstrated how an orbiting dense cloud 
of gas could also be detected through its time-dependent 
21 cm absorption. In this case we delineated a method 
to determine the spot’s orbital parameters from the ab¬ 
sorption spectrum. 

Our methodology does not require the system to be 
spatially resolved. While the angular size of Sgr A*’s 
photosphere is ~ 10 milliarcseconds, interstellar blurring 
could render even Sgr A* unresolveable. Interstellar ab¬ 
sorption would also deepen the absorption profile at low 
velocities. However, the absorption profile near the su- 
permassive black hole is wide, and at higher velocities 
(> 1000 km/s) the absorption should be purely from HI 
near the black hole. 

While we presented results for a homogeneous neutral 
disk with no limb darkening, our formalism is general 
and could be used for systems where Nh or I^o are non¬ 
trivial functions of the disk coordinates, or where the line 
function has a significant velocity width. 
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